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We study two aspects of the elasticity of smectic- A elastomers that make these materials genuinely 
and qualitatively different from conventional uniaxial rubbers. Under strain applied parallel to 
the layer normal, monodomain smectic- ^4 elastomers exhibit a drastic change in Young's modulus 
above a threshold strain value of about 3%, as has been measured in experiments by Nishikawa 
and Finkelmann [Macromol. Chem. Phys. 200, 312 (1999)]. Our theory predicts that such strains 
induce a transition to a smectic-C-like state and that it is this transition that causes the change in 
elastic modulus. We calculate the stress-strain behavior as well as the tilt of the smectic layers and 
the molecular orientation for strain along the layer normal, and we compare our findings with the 
experimental data. We also study the electroclinic effect in chiral smectic- A* elastomers. According 
to experiments by Lehmann et al. [Nature 410, 447 (2001)] and Kohler et al. [Applied Physics A 80, 
381 (2003)], this effect leads in smectic- A* elastomers to a giant or, respectively, at least very large 
lateral electrostriction. Incorporating polarization into our theory, we calculate the height change 
of smectic-A* elastomer films in response to a lateral external electric field, and we compare this 
result to the experimental findings. 

PACS numbers: 83.80.Va, 61.30.-v, 42.70.Df 



I. INTRODUCTION 



Smectic elastomers [l[ are rubbery materials that have 
the macroscopic symmetry properties of smectic liquid 
crystals @. They possess a plane-like, lamellar mod- 
ulation of density in one direction. In the smectic-A 
(SmA) phase, the Frank director n describing the aver- 
age orientation of constituent mesogens is parallel to the 
normal N of the smectic layers whereas in the smectic- 
C (SmC) phase, it has a component in the plane of 
the layers. Monodomain SmA elastomers are macro- 
scopically uniaxial elastomers, albeit with unusual me- 
chanical and electrical properties. For example, when 
strained along an axis normal to smectic layers, they ex- 
hibit a reorientation of smectic layers and an associated 
decrease in Young's modulus Yj| along this axis above 
a critical strain. This phenomenon, discovered experi- 
mentally by Nishikawa and Finkelmann (NF) Q, is the 
analog in smectic elastomers of the Helfrich-Hurault ef- 
fect in uncrosslinked smectics [1, Moreover, like ther- 
motropic chiral smectics, chiral Srm4 elastomers exhibit 
electrostriction, as depicted in Fig. [TJ in which smectic 
layer spacing decreases in response to an electric field in 
the plane of the layers. Early experiments @ produced 
a reduction in layer spacing by as much as 4% at field 
of 1.5 MV m _1 . More recent experiments Q produce a 
reduction of 1% in fields as high as 3 MV m _1 in films 
with one free edge to reduce mechanical stress. Even a 
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FIG. 1: Electroclinic effect in Sm^l* elastomers, (a) Without 
electric field, the elastomeric film has a thickness ho- (b) By 
application of a lateral electric field, a tilt angle <f)E is induced 
in a plane perpendicular to the field and the film thickness 
decreases to a value h. 

1% reduction in height is larger than that produced in 
most traditional actuators. 

In this paper, we analyze the experiments just dis- 
cussed using a phenomenological model of smectic elas- 
tomers we have recently developed p| Q. This model, 
which generalizes the standard Lagrangian approach to 
elasticity theory [l(| E|, brings together in a single rubric 
physics associated with the crosslinked network that de- 
fine the elastomer, the layering of the smectic phase, and 
the interaction between the nematic director and both 
the network and the smectic layers. It includes ener- 
gies that favor constant smectic layer spacing and that 
lead to the development of shear strain in response to 
stretch perpendicular to layers above a critical stress. It 
also includes symmetry-permitted interactions between 
c-ordering (molecular tilt relative to layer normal) and 
shear strain that cause shear to produce c-ordering and 
vice versa. Thus a strain along the smectic layer normal 
produces the c-ordering of a SmC elastomer, as well as 
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shear strain. We will analyze this instability of a SmA 
elastomer under strain parallel to the layer normal in 
terms of an induced transition to a SmC elastomer. Com- 
paring our numerical estimates for the stress-strain be- 
havior and the tilt of the smectic layers to the available 
experimental curves, we find convincing agreement. To 
analyze the electrostriction experiments, we add to our 
model the standard chiral couplings of the electric field 
to the Frank director and layer normal. The result is that 
an electric field in the plane of the layers produces both 
SmC ordering and shear strain linear in the electric field 
and a reduction in layer spacing quadratic in the field. 
Our numerical estimate of the magnitude of the reduc- 
tion in layer spacing is in accord with the experimental 
values. To our knowledge, experiments have not sought 
to detect the shear strain induced by the external electric 
field. 

Warner and Adams (AW) [12[ have recently developed 
a theory of the elastic properties of SmA elastomers that 
combines the neoclassical model of nematic elastomers 
1] with the compressional energetics of smectic layers. 
This theory assumes that the nematic director is locked 
to the layer normal and that the elastomer is incompress- 
ible. Our theory is phenomenological; it is based on an 
expansion, consistent with all symmetries, of the free en- 
ergy in powers of the Cauchy-Saint-Venant strain tensor 
[lOj l and the nematic director. The AW theory is semi- 
microscopic, and it applies, strictly speaking, solely to 
crosslinked networks of liquid crystal polymers. However, 
it provides a more realistic description of these particu- 
lar elastomers than does our phenomenolgical model if 
strains are large. We compare the AW theory and ours 
in some detail, and we explain how the former theory 
can be extended to allow for relative tilt between the 
director and the layer normal. Though the details of 
the AW approach and ours us differ considerably, except 
for the development of SmC order, they not surprisingly 
predict qualitatively identical results. A theory for chi- 
ral smectic elastomers was first set up by Terentjev and 
Warner [l3|, 03] using group theory. Our theory with po- 
larization generalizes the theory by Terentjev and Warner 
in a formalism that ensures invariance with respect to 
arbitrary rather than infinitesimal rotations of both the 
director and mass points. 

The outline of our paper is as follows: SectionHTlbrieflv 
reviews our model for smectic elastomers. Section IIIII 
treats the behavior of SmA elastomers under strain im- 
posed along the layer normal. Exploiting the model re- 
viewed in Sec. UH it derives predictions for the stress- 
strain behavior as well as the layer and the molecular 
tilt in response to strain. Section HVl incorporates polar- 
ization into our model of Sec. |TT] and considers the elec- 
troclinic effect and electrostriction. Section [V] presents 
concluding remarks. There is one appendix that com- 
pares our theory in some detail to the AW theory as it 
stands. Moreover, the appendix presents a generalization 
of the AW model that allows for SmC ordering. 



II. MODELLING SMECTIC ELASTOMERS 

Let us first review briefly our model for smectic elas- 
tomers. Smectic elastomers are, like any elastomer, per- 
manently crosslinked amorphous solids whose static elas- 
ticity is most easily described in Lagrangian coordinates 
in which x labels a mass point in the undeformed (ref- 
erence) material and R(x) = x + u(x), where u(x) 
is the displacement variable, labels the position of the 
mass point x in the deformed (target) material. La- 
grangian elastic energies are formulated in terms of the 
strain tensor u — ^(g — 6) where 6_ is the unit matrix, 

g = A T A is the metric tensor, and Ay = dRi/dxj 
are the components of the Cauchy deformation tensor 
A with i,j, k — x,y,z. The components of u are thus 
wy(x) = \ (Aj k A kj - Sy) = \ (diUj + djUi + diU k djU k ). 
Here and in the following the summation convention on 
repeated indices is understood. 

As in nematic elastomers, the Frank director n inter- 
acts with the elastic strain. The strain u is a reference- 
space tensor that transforms under the same reference- 
space operations that transform x. The director n, on 
the other hand, is a target-space vector that transforms 
under the same target-space operations that transform 
R. To create scalar invariants coupling u to n, we need 
to convert n to a reference-state vector n. This is ac- 
complished with the aid of the polar decomposition [l5| 
of the deformation tensor: A = Og 1 ^ 2 , where g 1 ^ 2 is 

a symmetric reference-space tensor and O = Ag -1 / 2 
is a rotation matrix. The components of O are Oy = 

Aifc[<7 -1 / 2 ]/cj where we use the notation that [M a ]y is the 
«j-component of the matrix M a for any square matrix M 
and exponent a (it is understood that My = [M]y). 

The left and right indices in Oy transform, respec- 
tively, under target- and reference-space operations, and 
O converts reference-space vectors to target-space vec- 
tors. Thus, we have n = On and n = T n. We will ex- 
press the vectors n and n in terms of their components c a 
and c a perpendicular to the initial anisotropy axis along 
z: rii = c a ,n z where n z = (1 — c 2 ) 1 / 2 w 1 — ^c 2 and 
similarly for hi, where we introduced the notation that 
letters a, b,c... at the beginning of the alphabet run over 
the directions x,y perpendicular to z. Reference-space 
vectors such as c and the vector e = (0, 0, 1) along the 
reference-space uniaxial direction [l6| can be contracted 
with the strain ity to produce scalars like c a u a bCb and 

We can now construct free-energy densities, which we 
will simply call energies, for the various contributions to 
the total energy / of a smectic elastomer. There is a 
contribution / net to / that is identical to the soft elastic 
energy of nematic elastomers including both the familiar 
terms quadratic in uy and c a . The characteristic ener- 
gies scales in this free energy are set by the volume bulk 
modulus B ~ 3 x 10 9 Pa, the network shear modulus 
/i ~ 10 6 , and moduli that scale as \i times a factor that 
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vanishes when the ratio p of the polymer step lengths 
parallel and perpendicular to the local director becomes 
unity. 

Smectic elastomers have a layer structure with a pre- 
ferred layer spacing whose magnitude depends on the 
angle between n and the layer normal N (or, equiva- 
lently between n and N). Smectic order is described by 
a complex mass-density-wave amplitude whose phase is 
<t> = Qo[Rz - U(R)\ = q [z + u z (x) - U(R{x))\, where 
U (R) is the displacement field of the smectic layers, and 
go = 2ir/d, where d is the preferred layer spacing when n 
is parallel to N. We consider only elastomers crosslinked 
in the smectic phase in which case, C/(R(x)) = u z (x), 
and <f> = qoz. The energy associated with changes in the 
smectic layer spacing is 

/layer = ^s m 9„ 4 [(n • V</>) 2 - q 2 } 2 , (2.1) 

where B sm is the smectic compression modulus with a 
value of order 10 7 Pa deep in the smectic phase though 
it vanishes as the nematic phase is approached. /i aye r 
favors a layer spacing of d' — d/ cos O, where as before 
is the angle between n and N. 

Finally, there are interactions favoring n parallel to N 
in the smectic- A phase and tilted relative to N in the 
smectic-C phase: 

/tut = \n sin 2 9 + X -v t sin 4 6 . (2.2) 

The energy /i aye r + /tilt is the generalization to smec- 
tic elastomers of the Chen-Lubensky model 17j for 



Srm4-SmC-nematic phase behavior in uncrosslinked liq- 
uid crystals. Its coefficient r t is linear in the deviation 
of the temperature T from the transition temperature T c 
between the SmA and SmC phases, r t = a\T — T c \. The 
coefficient Vt is essentially independent of temperature. 
Experiments by Brehmer, Zentel, Giefielmann, Germer 
and Zungenmaier (BZGGZ) [18] on smectic elastomers 
measured the values of the coefficients of a tilt energy 
similar to our f t nt (however, their tilt energy is somewhat 
more general in that it also accounts for polarization ef- 
fects which we will discuss later). A short calculation 
translates their experimental values into estimates for our 
a and v t : a ~ 1.3 x 10 4 Pa/K and D t ~5x 10 5 Pa. At 
a temperature of 20K above the transition temperature, 
r t is or order 10 5 Pa. More recently, similar experiments 
were performed on liquid smectics by Archer and Dierk- 
ing (AD) [l]|. Their results lead to a ~ 4 x 10 4 Pa/K 
and Vt ~ 10 6 Pa implying that r t ~ 10 6 Pa at 20K above 
the transition temperature. Though measured for liquid 
smectics, the AD values should be of some relevance to 
elastomeric smectics, and we will use them in the follow- 
ing in addition to the BZGGZ values in order to put the 
estimates that we are going to make on a broader basis. 

The energy of a smectic elastomer / = /net + /layer + /tilt 
can now be expressed as a sum of a harmonic uniaxial 
elastic contribution f unl depending only on Uy , a con- 



relevant nonlinear terms, a contribution f c describing the 
energy associated with the formation of non vanishing c- 
order described by c a , and a contribution / coup i arising 
from the coupling of Uij and c a : 



f /uni "b /nonlin "b fc "b /coupl- 



(2.3) 



The uniaxial elastic energy to harmonic order in strains 
is 



/uni — 9 ^1 



C 2 u z 



+ h C 3 u 2 



+ Ci u 2 ab + C 5 u z az 



where 



Uab = Uab~ h SabU c 



(2.4) 



(2.5) 



is the two-dimensional symmetric, traceless strain tensor. 
The nonlinear energy reads 



/nonlin = -B X U zz u\ z + B 2 (ul z ) 2 



The c-director energy is 

f c = lrcl + lv(cl) 2 , 

and the coupling energy is 

/coupi = Ai c 2 a u zz + A 2 c 2 a u it + A 3 c a u ab c b 
+ A4 c a u az + A 5 



(2.6) 



(2.7) 



(2. 



tribution / n 



depending only on Uy that collects the 



Table U reviews the relations between the original elastic 
constants of / net , /i ay or an d /tut an d the effective elastic 
constants featured in Eq. ([2^4]) and Eqs. ([23)1 to (JXSJ). As 
mentioned earlier, the parameter p appearing in Table U 
is the anisotropy ratio. In the work of NF, the sam- 
ples showed at the transition temperature a spontaneous 
stretch along the director of about 12%. When making 
estimates, we will view p ~ 1.1 as a typical value. Based 
on the relations given in Table U we deduce that the hi- 
erarchy of magnitudes of the effective elastic constants 
is: C 3 - 10 9 Pa, d ~ B 1 ~ B 2 ~ v ~ Ai ~ A 5 ~ 10 7 

Pa, C 2 ~ C 4 - C 5 - r - [i - 10 6 Pa, A 3 10 5 Pa 

and A2 ~ — 10 4 Pa. Our estimates for r and A4 depend 
noticeably on whether we use the BZGGZ or the AD val- 
ues: we expect r ~ 10 5 Pa and A4 ~ — 10 5 Pa based on 
BZGGZ and r ~ A 4 ~ 10 6 Pa based on AD. It is worth 
commenting on the sign of A4. Our estimates of its value 
have different signs for the two physical systems for which 
we have data. This is possible because the two contribu- 
tions to A4 have opposite signs. The second contribution, 
r t , which is positive, arises because O ~ c a + u az , and 
the I rt sin 2 O-term favors c a — —u az . The first term, 
~/i(p 2 — l)/p, describes the rotation of c a in response to 
the imposition of a shear u az . It is negative and it favors 
a c a with the same sign as u az . 

Several observations about / are in order. First, c a and 
u az are coupled at linear order via the A4 term. Thus, 
if a nonzero u az develops, it will be accompanied by the 
development of a nonzero c a and vice-versa. Negative 
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TABLE I: Contributions to coefficients in / from /net, /layer and ftut with a = /j,(p — l) 2 /(2p). 
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A4 favors mechanical tilt and mesogenic tilt in the same 
direction whereas positive A4 favors opposite tilt direc- 
tions. If, for example, u zz = and u az ^ 0, then c a 
relaxes to c a = — (^i/r)u az to minimize /. This leads to 
and effective or renormalizcd renormalizcd value of the 
shaer modulus C5, 
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If, to give another example, u zz = and c a 7^ 0, u az 
relaxes to u az — — X^/ {2Cc,)c a and r is renormalized to 
Tr = r — A|/(2Cs). Note that the value of the renor- 
malized elastic constants and does not depend on 
the sign of A4. Second, an externally imposed strain u zz 
reduces the coefficient of u 2 az via the B\ term. This leads 
to the elastomer version of the Helfrich-Hurault 0, [f| in- 
stability, which we will discuss in Sec. IIIII Finally, if c a 
becomes nonzero, there will be an accompanying decrease 
in u zz and the smectic layer spacing proportional to (? a 
via the c„u zz -term with positive coupling constant Ai. 
The reduction of layer spacing in response to an external 
electric field that tilts the director in chiral systems will 
be the subject of Sec. |IV] 



III. ELASTICITY OF SMECTIC-A 
ELASTOMERS UNDER STRAIN IMPOSED 
ALONG THE LAYER NORMAL 

Before we start with our analysis, it is useful to clar- 
ify the experimental conditions. In the experiments of 
NF [3j, the Srm4 sample is subjected to an external uni- 
axial stress normal to the smectic layers. This stress 
stretches the sample, producing a deformation tensor A 
with A zz > 1. If there is a transition to the SmC phase, 
the sample will undergo a shape change, as sketched 
in Fig. O from a rectangular parallelepiped to a non- 
rectangular one with two parallel non-rectangular faces. 
The external stress is applied along the space- fixed z-axis 
and causes the target space vector R to stretch along the 
z-axis. This defines a preferred orientation in the target 
space and therefore breaks rotational invariance in this 
space. In an idealized experiment, the external stress 
will stretch the sample parallel to the z-axis in the target 
space only, regardless of the shape of the non-rectangular 
faces. Thus, in the SmC phase, the configuration of the 
sample will be that (shown in Fig.^c)) of a simple shear 
in which A zx > and A xz = 0. In a real experiment, 
boundary conditions prevent a distortion with a spatially 
uniform value of A zx , and the system breaks up into op- 
positely tilted domains. We will ignore this effect in our 
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FIG. 2: Schematic representation of distortions in the xz- 
plane induced by the Sm^4-to-SmC transition: (a) undistorted 
SmA phase, (b) SmC phase with a symmetric deformation 
tensor A , and (c) SmC phase with A xz — and A zx > 0. 



analysis. Thus, the general form for the deformation ten- 
sor under uniaxial stress along the target-space z axis is 



A 



A = 



./•.;• 

A 







yy v 
A zx A 



(3.1) 



where A zx = in the SmA phase. This then implies a 
strain tensor 



1 / A L + Kx - 1 a m a„ 

1=0 A lv ~ 1 

1 V A ZX A ZZ A\ z - 1 



(3.2) 



Thus, the zz-component of the strain tensor, u zz — 
{A 2 ZZ — l)/2, depends only on the imposed strain defor- 
mation A zz in both phases. The formation of the SmC 
phase is signalled by the development of a nonzero shear 
strain u xz = A zx A zz /2. 

As we indicated earlier, the imposition of a sufficiently 
large positive strain u zz induces a transition to a new 
phase with a nonvanishing shear strain u az . Because of 
the linear coupling A 4 c u az between strain and c-director 
in jcoupi, director tilt develops along with shear strain, 
and the phase with nonzero u az can be viewed as a strain- 
induced SmC phase. Effective energies in terms of either 
Uaz or c a provide equally valid descriptions of the tran- 
sition. Here we will describe the transition in terms of 
an effective energy expressed in terms of c a . To arrive 
at this energy, we must relax the variables Uu, u a b and 
u az to their equilibrium values in the presence of nonva- 
nishing c a and u zz . The equations for df/duu = and 
df /du a b — are linear and are easily solved for Ua and 
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Uab- 



Uab 



C2 _ ~2 
^ u zz Tc ai 
^3 

-L) (c a C b - \5ab^ c ) , 



where 



t = A2/C3 and u = A 3 /(2 C 4 ). 
The equation of state for u az is nonlinear: 
Of 



(3.3a) 
(3.3b) 

(3.4) 



du ai 
where 



= 2C 5 (u zz )u a b + AB 2 u bz u az + A 4 (m Z2 )c q , (3.5) 



C 5 (u zz ) = C 5 - Biu zz , 
A4 («zz) = A 4 + X 5 u zz . 



(3.6a) 
(3.6b) 



Since we are interested in properties near the transition, 
we solve Eq. (|3.5|) for u az to third-order in c a : 



where 



p{u zz ) = 



b'-ai 



2B2-— — -cc, 

C 5 (u zz ) 

Xi(u zz ) 



2c b (u zz y 



(3.7) 



(3.8) 



Inserting Eqs. (|3.3a[) . (|3.3b[) and (|3.7p into the free energy 
of Eq. H'2.3p . we obtain an effective energy expressed as a 
function of c a and u zz (the latter being fix externally), 

Us = |Cf t4 + |r fl (u„)^ + i Ur ( u ^)c2c2, (3.9) 
where 

Cf = Ci - C 2 2 /C 3 , (3.10a) 

rn(u zz ) = r + 2(Ai - C 2 r)u zz - 2C 5 (u zz )p 2 (u zz ) , 

(3.10b) 

= « - 2C 4 w 2 - 2C 3 r 2 + 4B 2j0 4 ( Uz2 ) . (3.10c) 

We can now use f e g to analyze the transition to the 
SmC phase induced by the uniaxial strain u zz . The term 
2C 5 (u zz )p 2 (u zz ) = Xl(u zz )/2C 5 (u zz ) in r R (u zz ) diverges 
at u zz = C5/B1 and guarantees that rn(u zz ) will pass 
through zero upon increasing u zz . u c zz , the critical value 
of u zz at which = 0, is given by 



= A" 



S - a/S 2 - 2rCf A 



(3.11) 



with as given in Eq. (|23|) . and where A = 4i?i (rC 2 - 
Ai) - A^ and S = rB x + 2tC 2 C 5 - 2C 5 \i + A 4 A 5 . Both, 
the BZGGZ and the AD values lead to u c zz w 0.03 in 
good agreement with the experimental estimate of 0.034 
[3j. In the vicinity of the critical point determined by 



u c zz , we can expand rn{u zz ) and gn(u zz ) to lowest order 
in 



5u zz = u zz - u c zz : 



tr{u zz ) = -b5u zz , (3.12a) 
vr{u zz ) = v R (u c zz ) = v R , (3.12b) 

where b = —drR(u zz )/du zz \ u ^ =u c z > 0. Choosing our 
coordinate system so that the c-director induced by the 
u zz strain lies along the rr-direction we obtain = and 



c° 



for u zz < u zz , 

±s^/Su zz for u zz > u c zz , 



(3.13) 



where s = yJbjv c R . The director in the SmC phase 
under a symmetric shear, i.e. in the geometry shown 
in Fig. H)Jb), is therefore n = (sin a, 0, cos a), where 
sin a = (P x . Using the BZGGZ and the AD values, we 
estimate s w 0.75 and s « 0.8, respectively. 

From the equilibrium values of the c-director, it is 
a straightforward exercise to determine the equilibrium 
strain tensor vP from Eqs. (|3.3ap . (|3.3b[) . and (|3.7[) . Con- 
sistent with Eq. (|3.2p . there is no shear induced by u zz 
in the icy-plane, u® y = 0, and there is no shear induced 
in the yz-plane, vP yz = 0. There is, however, a shear in 
the xz-plane, which according to Eq. (|3.7[) and (|3.13[) is 



for u zz < u c zz , 
for u zz > ut, , 



(3.14) 



to lowest order in 8u zz with p c = p(u zz ). The BZGGZ 
and AD values, respectively, lead to p c ~ 1.2 and p c m 
1.1. 

The equilibrium values for extensional strains follow 
from Eqs. (|3.3a[) . (|3.3b[) and Ua = u aa + u zz . To lowest 
order in Su zz , 



9i 

C2 

c 3 



for u zz < ul 



9i 
c 3 



+ (t + Lu)a 2 



5u z 



for u zz < ul 



(3.15a) 



(3.15b) 



and u yy is given by Eq. (|3 . 1 5[) with uj replaced by —to. 
The parameter uj = \z/{2C§) = —(p,(p— 1) / p)[r t + p(p — 
l) 2 /p] is negative and of order —0.05. Figure [3] depicts 
the equilibrium strain components on u zz schematically. 

The quantities c° and u® z determine the equilibirum 
angle 0° between the layer normal and the director. To 
lowest order, 



c° 



Pc]c° x . 



(3.16) 



The coefficient 1 — p c of (P x on the right hand side is 
negative and of the order of —0.2 or —0.1, depending on 
whether we use the BZGGZ or the AD values. The small 
but non-zero 1 — p c implies that the tilt angles of the layer 
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FIG. 3: Dependence (schematic) of the non-zero equilibrium 
strain components on u zz . 



provided that Cf- ^ D\ , and hence the Young's modulus 
Yy is given by 

v _ / Ci for w zz < u c zz , , . 

in this case. Adams and Warner in fact find Cf 1 = Z?i 
in their model, and we obtain the same result in our 
treatment if we lock n to N and set / = /net + /layer- This 
limit does not follow cleanly from the model in Eq. (|2.3|) 
and the parameters of Table |l] because that model does 
not include all of the fourth-order terms in the free energy 
of the original model defined / = / nct + /i aye r + /tilt ■ If 
Cf = Di, then cr™ g = C^u c zz A zz so that the leading 
terms for u c zz -C Su zz <C 1 are 



normal and the director are distinct and that, however, 
their difference will be small for u zz not too far above 
u zz . We will revisit the tilt angles of the layer normal 
and the director further below. 

Now we turn to the effective elastic energy density as 
a function of u zz only when the other strain components 
and c have relaxed to their equilibrium values. Plugging 
c°, into Eq. (|3.9j) . we obtain to quadratic order in u zz 



for 



< u. 



(3.17a) 



and 



/ = ±Cf {u c zz + 5u zz f - \D X (Su zz ) 2 for u zz > u c zz , 

(3.17b) 

where D\ = (1 / 2)b 2 / v c R . Note that the two branches 
(|3.17a|) and (|3.17b|) match at U zz as they should. 

Next, we address the relation between stress and 
strain. The engineering or first Piola-Kirchhoff stress 
measured in experiments is c™ 6 = A^cry, where <7kj 
iif /dukj is the symmetric second Piola-Kirchhoff stress 
tensor. The system reaches equilibrium at fixed u zz with 
respect to all other components of My, i.e., = for 
ij =/= zz. Thus, the zz component of the engineering 
stress tensor is A zz cr zz , and 



CfA, 



for u zz < u c . 



(3.18a) 



and 



o c zz & = C?u c zz A zz + [C* - D 1 ]A zz Su zz for 
For the deformation tensor of Eq. (|3.ip and small strain 



A 



ZZ 

u° zz is small as our estimates indicate, 
a™ s w Cfu zz for u zz < u c z 
For u zz > u c zz and small 5u zz , 



1 



it. 



(3.18b) 



Thus, if the critical strain 



(3.19a) 



o£* « Cf u zz + [Cf - D t ]5- 



u. 



for u zz > u% 



(3.19b) 



for 



implying that 



for 
for 



Z2 < 
2Z > < 



(3.21) 



(3.22) 



if Cf = D x . Equations (j3TT9"|) and (j3~^T|) reveal that cr™e 
as a function of u zz will consist of two straight lines of 
different slope (Young's modulus) meeting at u zz — u c zz . 
If Ci = D\, then the ratio of the Young's modulus for 
u zz > u c zz to that for u zz < u zz is simply u c zz . Our full 
equations (|3.18|) show that there are nonlinear corrections 
to cr™ s both above and below u c zz . Since the experimental 
value of u zz ~ 0.034 is much less than one, the nonlin- 
ear corrections below u zz are unimportant. Those above 
u zz vanish or are unimportant if C-p and D\ are equal or 
nearly so. Figure [4] shows the linearized stress-strain be- 
havior along with experimental data. Both, Eq. (|3. 19b[) 
and (|3.2ip are used in the plots and both produce excel- 
lent agreement with the experimental data. 

For completeness we note that, if we use our strain- 
only model for the Srm4-to-SmC transition Q, we find, 
to leading order in 5u zz , essentially the same equilibrium 
strain tensor, the same effective elastic energy, and the 
same stress-strain relation as AW. Obvious differences 
reside in the different definitions of the elastic constants 
used in both models. These differences are of course qual- 
itatively unimportant. 

Next we calculate the tilt angle ip of the layer normal 
in response to u zz . As illustrated in Fig.O the angle that 
N makes with the z axis is 



f = + 1, 



(3.23) 



where j3 is the angle that N makes with the z-axis and 
7 is the clockwise angle about the y-axis through which 
the sample under symmetric shear has to be rotated to 
bring it into the configuration with A xz = 0. In terms 
of the components of the symmetric deformation tensor 
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0.05 0.1 0.15 0.2 



FIG. 4: The engineering stress <r™ g as a function of strain 
u Z z of a monodomain SmA elastomer. The squares symbolize 
experimental data by Nishikawa and Finkelmann 3]. Above 
a threshold strain of about 3% the elastomer is in a SmC- 
like phase with an elastic modulus much lower than that of 
the original Srm4 phase below the threshold. The continuous 
lines stem from Eqs. (I3.19a|) and (|3.19b|) with u c zz = 0.034, 
Cf = 3.5 x 10 7 Pa, and D x = 3.223 x 10 7 Pa. The dashed 
lines stem from Eqs. (|3.19a[) and (pOT]) with u% z = 0.0343 
and Cx — 3.5 x 10 7 Pa. The continuous and dashed lines 
lie almost on top of each other and are therefore hard to 
distinguish. 
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FIG. 5: Dependence of the tilt angle tp on u zz . The squares 
symbolize data taken from Ref. [3j. The continuous curve 
corresponds to our analytic result (|3.25p . 

Fitting our result (|3.25|) to the experimental curve, we 
obtain u c zz = 0.034, A x = 90 and A 3 = -49. The fitted 
analytical and the experimental curve agree remarkably 
well, and the best-fit value A\ = 90 is certainly of the 
same order of magnitude as our estimate A\ 100 

Finally, we calculate the tilt angle p n of the mesogens 
in response to u zz . The angle that the target-space di- 
rector makes with the z axis is, as depicted in Fig. [3J 



A g , these angles are given by 



j3 = tan" 



7 = tan 



/ A 

-1 / lv Szx 



A° Sx J 1-|(1 + <7 a /C 3 )«S 



Szx 



A 



Szz 



i + «S 



(3.24a) 
(3.24b) 



We find that tp vanishes as 8u zz — > 0, and near 5u zz = 0, 
its expansion to order (Su zz ) 3 ^ 2 is 

_ f for u zz < u c zz , 

tf ~\A 1 {5u zz ) 1 ' 2 + A 3 {8u zz f' 2 for u zz > u\ z . 

(3.25) 

The value of A\ is determined entirely by the lowest order 
expressions for v? xz [Eq. !pTH]l ] and u° xx [Eq. (f3~15)) ]. 



At = sp c 



2 + \{l-C 2 /C 3 )u zz 
[l+u zz ][l-Ul + C 2 /C 3 )u zz ] 



(3.26) 



A 3 has contributions both from the expansions of tan -1 
to third order in its arguments, and from (5u zz ) 3 / 2 con- 
tributions to u xz and u xx , which we have not calculated. 
Both, the BZGGZ and AD values, lead to At w 1.8. 
Next, we compare our result (|3.25p and our estimate for 
A\ with the experimental data by NF, see Fig. For 
this comparison it is important to note that the expan- 
sion (|3.25[) with A\ given by Eq. (|3.26[) expresses ip in its 
natural unit, radian. The experimental curve by NF as 
shown Fig. ([5]), on the other hand, measures p in degrees. 
Converting the above estimate for A\ from radian to de- 
grees we get the estimate A\ w 1.8 x 360/(27r) w 100. 



■a- j. 



(3.27) 



where a = sin^ 1 c° 
the vicinity 
is 



in ' c x . Like tp, tp n vanishes as Su zz — > 0. In 
of Su zz — 0, its expansion to order (Su zz ) 3 ^ 2 



for u zz < u c zz , 



^ 1 A nA {5u zz fl 2 + A n , 3 (5u zz f/ 2 for u zz > u% z , 

(3.28) 

with the value of ^4 n .i being determined entirely by 
the lowest order expressions for c x [Eq. (|3.13[) ] and u xz 
[Eq. ([331], 



A T 



Spc 



1 



(3.29) 



^■n,3 has contributions both from the expansions of sin -1 
and tan" 1 to third order in their arguments, and from 
(5u zz ) 3 / 2 contributions to c x and u xz , which we have not 
calculated. From the BZGGZ and AD values we obtain, 
respectively, approximately A n i w 85 and A n i w 90 
when tp a is measured in degrees. 

In their experiments, NF did not seek to produce 
curves of tp a versus u zz . From measuring the layer thick- 
ness as a function of u zz and by interpreting their X- 
ray reflection patterns, however, they conclude that the 
layer normal and the director remain parallel above the 
threshold strain and that the reduction of Yj| above the 
threshold stems form a strain induced breakdown of the 
smectic layers into pieces with local Srm4 order. Our the- 
ory, on the other hand, predicts that the angles tp and 
tp a are different but of the same order of magnitude, and 
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that the reduction of Yj| results from shear and SmC-like 
ordering produced by sample extension. If the difference 
between (p and tp a is small, it can be difficult to distin- 
guish between SmA and SmC order experimentally. If 
this is the case in the experiments of NF, then the dif- 
ferences between the experimental findings and our the- 
oretical predictions for the tilt angles may lie within the 
experimental error. 

Above, we have at several occasions briefly commented 
on the relation of our work to that of AW. A more de- 
tailed comparison of the two approaches can be found in 
the Appendix. 



IV. ELECTROCLINIC EFFECT IN SMA* 
ELASTOMERS 

In a chiral smectic elastomer an external electric field 
E can couple to both the elastic and the orientational 
degrees of freedom. First, we consider the former cou- 
pling. E, like the director n and the layer normal N, is 
a vector in the target space and hence cannot be cou- 
pled directly to the reference space tensor u. By now, of 
course, we know how to overcome this problem: we can 
use the polar decomposition theorem to switch from E to 
its reference space counterpart E via the transformation 
E = O t E and then construct form E and u couplings 
that are manifestly invariant under simultaneous rota- 
tions of the sample and the field in the reference space. 
Since we are interested in SmA* elastomers, the couplings 
have to conform with uniaxial symmetry in the reference 
space which leads us to 



Iea = -«i E a u ab E b 



v 2 u zz E z , 



(4.1) 



where v\ and i>2 are coupling constants. In Eq. (|4.ip 
we have omitted a term of the type u aa E 2 to keep our 
discussion as simple as possible. This omission does not 
affect our results qualitatively. 

Now we come to the coupling between the electric field 
and the orientational degrees of freedom. In a chiral 
smectic a spontaneous polarization P occurs if the di- 
rector n and the layer normal N are neither parallel nor 
perpendicular: 



n x N(n • N) . 



(4.2) 



Thus, in the presence of E, the elastic energy density of 
the elastomer has an extra contribution 



/ B)2 ~ E • n x N (n • N) . 



(4.3) 



In order to combine ]e.2 with the remaining terms of 
the elastic energy density we recast Je,i in the reference 
space to obtain 



f Ef2 = K E-nx N(n-N) 



(4.4) 



where K is a coupling constant. The experimental find- 
ings by BZGGZ imply K ~ 10" 3 Pa(m/V), whereas the 



value of K implied in the results by AD is one order of 
magnitude smaller, K ~ 10~ 4 Pa(m/V). 

In the typical experimental setup, the normal of the 
elastomeric film is along the z direction and the elec- 
tric field E is along the y direction. If we assume that 
the experiments correspond to the geometry shown in 
Fig. [H then A° a = but A xz ^ 0, and the only nonzero 
components of A are are A® x , A yy , A zz , and A xz , im- 
plying that that O xy = O yx = O yz = O zy = 0. Thus 
Ei = OfjEj = Ei, and E = (0,^,0). In addition, 

Ni = [g-^hig- 1 ]^ 2 « (-««*, 1- With these 

expressions, the above coupling energy densities simplify 
to 



Je,i = —Vi u yy E 2 , 
f E ,2 = ~KE(d x + 



(4.5a) 
(4.5b) 



In Eq. (|4.5bp we have dropped higher order terms in c a 
which are inconsequential for the leading behavior since 
we are interested in the SmA* and not in the SmC* 
phase. For the same reason, we need, when we com- 
bine Je,i and ]e,2 with the remaining parts of /, to keep 
in f c and / COU pi only terms up to second order in c a . To 
keep our discussion as simple as possible, we will focus 
in the following on the incompressible limit. If strains 
are small, as we assume, we can implement this limit by 
setting uu = 0. In accordance with the given geometry, 
we set u xy = Uy Z — c y — 0. Collecting all parts of the 
elastic energy density we then obtain 

/ = Ami + /nonlin + fc + /coupl + JE,l + Je,2 

= 5 (Ci + C 4 ) u 2 zz + 2C 4 (u 2 xx + u xx u zz ) + C 5 u 2 az 



+ \rc 2 x + |(2Ai +A 3 )c 



B 1 u zz ul, 



+ A3 c 2 U xx + A4 c x u xz + \^u zz u xz c x 
- K E (c x + u xz ) + vi {u xx + u zz ) E 2 



(4.6) 



where we have eliminated u yy by using ua = and where 
we dropped the B 2 (u 2 z ) 2 term. 

Next, we determine the equilibrium values of the strain 
components and the c-director in the presence of the elec- 
tric field. The equilibrium values c x and u° zx , which are 
at leading order linear in E, are 



2C 5 r R 
A4 — r 
~2C 5 r R 



KE 



r- A 



2C5 — A4 



4 ~ c 



(4.7a) 
(4.7b) 



where tr = r — X 2 /(2C^). The diagonal components of 
vf-j can also be calculated by minimizing the free energy 
of Eq. (14. 6|) . The results are that u° zz and u® x both have 
contributions proportional directly to E 2 and contribu- 
tions proportional to (c x ) 2 . In addition u° zz has a contri- 
bution proportional to —B{u xz ) 2 . All of these terms are 
proportional to E 2 . 

The relative height change 



Ah/ h = (h-ho)/h 



(4.8) 
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of the film is 



1 

1 



&h/h = -^[Xi(c a x f + \ 5 c° x u xz - B 1 (u° xz f] - ^E 2 



2Ci 



[{%+u xz y-4(u xz y 



2d' 



E z 



(4.9) 



where we used B sm 3> /i, r t to obtain the last form (which 
sets /layer = 0). The electrostriction coefficient is defined 
via Ah/ho = —a c \E 2 . Inserting Eq. (|4.7[) into Eq. 
we obtain 



Vl 



2(2p-l) f2C 5 -X 4 

A 



2d (p+1) 2 V 2C 5 r R 
i 2p - 1 K 2 



2Ci 2 p 2 



(4.10) 



Finally, we compare our theoretical findings to the ex- 
perimental results. In their experiments, Lchmann et al. 
applied alternating lateral voltage, U ac (t) — U ac cos(ujt), 
so that E = E ac (t) = U ac (t)/d, where d is the width of 
the sample, and they measured the first and the second 
harmonic of Ah. Rewriting Eq. (14. 9f> as 



Ah= hoa A c \^ c [l + cos(2^)] 
4d 2 



(4.11) 



our model predicts that the amplitude Ah 2 of the second 
harmonic is 



Aho 



hp a c i Ul 
Ad 2 



(4.12) 



In Fig. [5] our result (|4.12[) is compared to the exper- 
imental data of Lehmann et al. Taking the values 
ho = 75 ± 5 nm and d = 1 mm of Ref. [6| and setting 
del = 7.147 x 10~ 14 (m/V) 2 , we obtain an excellent fit 
between our theory and the experimental curve for the 
second harmonic. Likewise, we can fit our theory easily 
to the experimental curves of Koler et al. 0]. We re- 
frain here to show the corresponding curves in order to 
save space. Beyond producing fitting curves, we can ex- 
tract from our theory an estimate for the value of a e \. 
Let us start with the first term in the second line of 
Eq. (|4.10p . The coefficient vi should be of order e Xo 
where e = 8.85 x 1CP 12 Pa (m/V) 2 is the permeability 
of the vacuum and xo ~ 2.5 is the electric susceptibil- 
ity. With these values, the term containing vi is of or- 
der 2 x lCr n /4 x 10 7 w 0.5 x 10~ 18 (m/V) 2 , which is 
orders of magnitude smaller than the experimental val- 
ues. Now, we turn to the second term in the second 
line of Eq. (|4.10p . In accord with this term, the ex- 
periments of Lehmann et al. and Koler et al. found the 
largest values of a D \ for temperatures near T c . Based on 
the BZGGZ and the AD values, we estimate r t to be of 
the order r t ~ 10 4 Pa and r t ~ 10 5 Pa, respectively, 
for T within a range of a few degrees about T c . Re- 
calling the experimental values for K and that p w 1, 
were are led to a e i « 5 x 10~ 15 (m/V) 2 for the BZGGZ 
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FIG. 6: Hight change of a Sm^4* elastomer in response to 
a alternating lateral voltage U ac (t) — U ac cos(ojt) as a func- 
tion of the amplitude U ac . The data points are taken from 
Ref. 0], where a sample of thickness ho — 75 ± 5nm and 
width d — 1 mm was used and the frequency of the voltage 
was uj = 133 Hz. There is a small piezoelectric contribution, 
i.e., a contribution that is linear in [/ ac or respectively the 
amplitude of the electric field Ea. c (t) = Ua. c {t)/d. In contrast, 
the contribution proportional to the square of the field, i.e., 
the effect of electrostriction, at the second harmonic is large. 
The continuous curve stems from our theoretical result (|4.12p 
with del = 7.147 x 10~ 14 (m/V) 2 . 



values and a e \ w xl0~ 18 (m/V) 2 for the AD values. Per- 
haps not surprisingly, the latter estimate, being based 
on experimental data for liquid smectics, turns out poor. 
The estimate a c \ ~ 5 x 10~ 15 (m/V) 2 is in much bet- 
ter agreement with the electrostriction experiments. It 
is about one order of magnitude smaller than the value 
a e i = 7.147 x 10~ 14 (m/V) 2 obtained by fitting the data 
of Lehmann et al. and it is of the same order of magni- 
tude as the value a \ = (1 ± 0.2) x 10~ 15 (m/V) 2 quoted 
in Ref. 0. 



V. CONCLUDING REMARKS 

We have investigated the behavior of monodomain 
Smj4 elastomers under strains u zz along the smectic layer 
normal. When exceeding a threshold value u c zz , these 
strains induce a transition to a SmC-like state. This ef- 
fect is accompanied by a tilt of the mesogenic component 
that sets in at the same threshold. Due to the transition, 
the Young's modulus Yji is qualitatively different from 
that of conventional uniaxial elastomers: the slope of the 
stress-strain curve changes from high constant value for 
u zz < u c zz to a much lower constant value for u zz < u zz . 
Our results for the stress-stain behavior as well as for 
the layer tilt agree nicely with the available experimental 
data. 

Moreover, we studied the electroclinic effect in Srm4* 
elastomers in an external, lateral electric field. Here, it is 
the external field that induces a transition to a SmC*-like 
state. Since this state is sheared in the plane perpendic- 
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ular to the field, the height of a sample is lower than 
in the initial Srm4* phase. Our theory finds, in abso- 
lute agreement with the experimental evidence, that his 
height change is proportional to the square of magnitude 
of the external field, i.e., that the mechanism at work 
is the so-called electrostriction. Our numerical estimate 
for the electrostriction coefficient is in accord with the 
available experimental values. 
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NOTE ON THE WORK OF ADAMS AND 
WARNER 

Here we will discuss the connection between the AW 
theory and our theory for SmA's in some more detail to 
understand the agreement in the results more systemat- 
ically. Moreover, we will generalize the AW theory to 
allow for relative till between the director and the layer 
normal. 



1. Comparison to the work of Adams and Warner 

AW extended the neoclassic model of rubber elasticity 
by formulating an elastic energy density 



/ 



AW 



i/xTr 



44 4 T I _1 



\B{d/d Q - 1) 



(i) 



The first term on the right hand side of Eq. (fTJ) is the 
usual trace formula of the neoclassic model, [i is the 
shear modulus. 



6+(p- l)n n 



(2) 



is the so-called shape tensor describing the conformations 
of polymeric chains before deformation and 



(1 — p 1 )nn 



(3) 



is its counterpart after a deformation has been applied. 
As in the main text, p denotes the anisotropy ratio of 
the uniaxial state, which we assume to be prolate, p > 
1. The AW model assumes incompressibility, i.e., the 
deformation tensor is subject to the constraint det A = 1. 
The new element in the theory of AW, the second term 
on the right hand side of Eq. |TJ, describes changes in 
the spacing of smectic layers. B is a layer compression 
modulus, do and d are respectively the layer spacing 
before and after deformation which are related via 



d 
d 



1 



(A- i ) T n 



(4) 



In the following, we chose n = (0,0,1), and we set 
n = (sin#, 0, cosd). Assuming a deformation tensor of 
the form shown in Eq. (|3.ip . the elastic energy density 
becomes 

/4m{aL + aL + a-a- + ,-aL 
(i-p^XML-aL + aL)^ 

- 2(1 - p~ 1 )A xx A zx sin 9 cos 9} 

2> 



, I A ZZ A XX 
o — — — 1 



A 2 



(5) 



where A yy has been eliminated via the incompressibility 
constraint, and where b = B/fi. For studying the re- 
sponse of an SmA elastomer to an imposed deformation 
A zz , AW assume that n and N are locked, which amounts 
to setting 



A, 



A, 



vAL + Ajc v / Axk + Afa 
Then, the AW elastic energy density becomes 



(6) 



/aw - tja4 AL + A2 A2 + A 2 za 



^-xx-^zz 



(AL+pApA 2 
A 2 +A 2 

XX 1 zx 



V ^lx + 



(7) 



In order to make contact to our Lagrangian theories, 
we re-express the components of the deformation tensor 
in terms of the components of the strain tensor. From 
the definition of the strain tensor and (|3.1[) it is straight- 
forward to see that the respective components are related 
by 



Ay 

A, 
A, 



\Jl + 2llyy , 

VTT2uT z , 

^xz I A zz , 



A xx = + 2u xx 



A? 



(8a) 
(8b) 
(8c) 

(8d) 



Inserting Eq. ([8]) into the elastic energy density ([7]) to 
eliminate the A^- and then expanding in powers of the 



Uij we obtain 



/aw = A* 

4 + 6 



4 + 6 



P .2 



2 ■ 2 u zz + 2u zz u xx + 2u xx + -u- zx 

4 + 6 

Uz Z u 2 zx - (1 +p)u xx u 2 zx + — — u zx \ , (9) 



where we have dropped higher order terms which are in- 
consequential for the argument here. 

Next, we revisit elastic energy density / in Eq. (|2.3p . 
To obtain a model in terms of strain only, we integrate 
c a out of /. This procedure leads to 

/ = Ami + Di U zz U 2 az + D 2 UuU 2 az + D 3 UabUazUbz 

+ G(u 2 az ) 2 , (10) 
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with G5 renormalized to G,f [see Eq. (I2.9|) ] in / un i, with 
D m = A m A 2 /r for m = 1, 2, 3, and G = ( 5 /4 + r/2)A^/r 4 . 
Then, we carry out the following steps: (i) we implement 
the incompressible limit via the constraint uu = so that 
the terms featuring G2, G3, and D2 are absent, (ii) we 
eliminate u vy with help of uu = 0, and (iii) we assume a 
deformation tensor of the form (|3.ip so that u xy = u yz = 
0. Then, / reduces to 



2C 4 u z 



+ 2dui x + a 

it 2 + Gui^ 



(11) 



which has, besides lacking the inconsequential constant 
term, exactly the same form as /aw- Comparison of the 
individual terms in Eqs. ^ and (fTTj) reveals that the 
elastic constants in the 2 models are related by C\ = 
M(3 + b),C 4 = Cf = pn, D x = M (p - 3 - 6)/2, D 3 - 
and G = /<t(4 + b)/8. Note that AW work with 
the values 6 = 5 and p — 2 which corresponds to negative 
values of D\ and D3. 

The effective elastic energy density in terms of u zz only 
produced by Eq. (fTTj) is essentially the same as the in- 
compressible limit of Eq. (|3 . 1 7[) . Thus, in the end, the 
AW model and our models lead, at least to leading order, 
to the same stress-strain curve, cf. Fig. @] 



2. Extending the model of Adams and Warner 

Now we generalize the AW model by avoiding the as- 
sumption that the layer normal and the director are 
locked. The full elastic energy density (JSJ) rather than 
the AW elastic energy density ([7]) will be the vantage 
point for the following considerations. In the remainder, 
we abbreviate A = A zz for notational simplicity. 

We can reduce Eq. J5| to a function of A, A xx , and 
A zx only by minimizing over 9. The result is 



sinf 



cos 9 



1 - 



B 2 
l + A 2 



62 

A 2 



where 



2 

XX 



A = pA 2 + A 2 ZX - A; 
B = 2A XX A ZX . 

The free energy minimized over 9 is then 

fm = 1^(1 + P ) A 2 + (1 +p- 1 )A 2 xx + 2 A 1 A 
+ (1+ P - 1 )A 2 ZX 



(12a) 
(12b) 



(13a) 
(13b) 



2 

XX 



(1 - p- 1 ) ^(rA 2 A 2 XX + A 2 ZX ) 2 + 4ALAL 



(14) 



To determine instabilities toward the development of A z 
we expand f m to second order in A zx : 



fn 



A 



2 

XX 



1 / A 2 

2 f \pA 2 -A 2 x 



2l 

XX \ 



A 2 



A- 2 A 
A(A 



b(A-l) 2 ] 



2b- 



A 2 

XX 



AL 



(15) 

= 0. 
1. 



Now we minimize this function over A xx at A zx -- 
The result is A xx = A -1 / 2 . Thus, when A = 1, A xx 
A 2 X — A 2 , and the coefficient of A zx is zero, i.e., there 
is no restoring force for infinitesimal shear strains in this 
model of smectics. The origin of this behavior is the 
invariance of both terms in the free energy with respect 
to rotations. Semi-soft terms need to be added to endow 
the smectic layers with a preferred direction relative to 
the crosslinked matrix. Next, we consider what happens 
when A > 1. The coefficient of A 2 X becomes 



A, 



1 , . A 2 + A 



1 



(16) 



This terms is negative for all A > 1 because 6^1. Thus, 
in the absence of a semi-soft term, there is no threshold 
for the production of shear strain in response to an im- 
posed strain along the z direction. 

As we have just seen, there is an is an instability to 
transfer shear for all A > 1 in the original AW theory 
if the nematic director is allowed to relax. Thus, the 
model elastic energy density ([5]) is insufficient to produce 
stress-strain curves as measured by NF. The problem is 
that this model is invariant with respect to simultane- 
ous rotations of the smectic layers, the nematic director, 
and A in the target space. To break this invariance, we 
introduce the semisoft energy 

/semi = ^oTr[(| - n n )A T nnA] (17) 
= -,ua(A xx sin0 + A zx cos6) 2 . (18) 
The complete energy is then 

h = f + /semi- (19) 

Expanding to second order in 9, we obtain 

/ = ^m|a 2 + A 2 XX + A,; ! A 2 + {p- 1 + a)A 2 zx 
- 2(1 - p- 1 - a)A xx A zx 9 

+ [(1 - p-^ipA 2 - A 2 , + A 2 J + a(A 2 xx - A 2 ZX )}9 2 

2^ 



AA T 



1 



Minimizing over 9, we obtain 
a {1 - p- 1 - a)A xx A z 



(1 - P -i)(M 2 - AL + AL) + «(AL - AL-) 



(20) 



(21) 
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and a relaxed free energy to harmonic order in A zx of 
1 



-a[ P A 2 + A 2 XX + A- 2 A, 2 } + b(A - 1)< 



1 

2^ 



(p 1)(A 2 - AL) + M(P 1)A 2 + AL) 



(p-i)(M 2 -AL) + «al 



2& 



A(A-l) 



AL 



A 



(22) 



The term of order in A zx has the same form as it did 
in the absence of the semi-soft term. Thus, as before, 
A xx = A -1 / 2 , and the coefficient of A 2 ZX is 



(p - 1)(A 3 - 1) + ap [(p - 1)A 3 + 1] 
(p- l)(pA 3 -i) + ap 



- 26A 2 (A - 1) 
If A = 1 + u zz , then to lowest order in u ZZl 

„2 



(23) 



A = 



ap 



(p — l) 2 + ap 



2bu 



3(p-l)(p(a-l) + l) 2 
({p-l) 2 +paf 



(24) 



and thus there is a nonzero the critical value of the strain 

u zz , 



p 2 a [(1 — pj 2 



pa\ 



zz 26 [(1 - pY + pa] - 3(p - l)[p(Q - 1) + l] 2 ' 

(25) 

in accord with the experiments of NF and our Lagrangian 
theory. Note that this critical strain is linearly propor- 
tional to a at small a and vanishes as a — > in agreement 
with our previous observation that the elastomer is un- 
stable to shear at infinitesimal strain along z if it is soft 
rather than semi-soft. Having u zz , we can proceed to 
calculate the entire stress-strain curve. The result is in 
full qualitative agreement with the experimental curves 
and the predictions of our Lagrangian theory. To save 
space, we leave the remaining steps as an exercise to the 
reader. 
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